program define mynormal_rdu0
	version 9.1
	args todo b lnf
	tempvar mu lnsigma sigma
	mleval `mu' = `b', eq(1)
	mleval `lnsigma' = `b', eq(2)
	quietly {
		gen double `sigma' = exp(`lnsigma')
		replace `lnf' = ln( normalden($ML_y1,`mu',`sigma') )
	}
end
